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Abstract. For given matrix A g 2 dxn , the set Pj, = {z : Az = b, z £ Z^} describes the 
preimage or fiber of b £ Z d under the Z-linear map fa : Z™ — » Z d , x i— > Ax. The fiber p, is 
called atomic, if Pj, = P^ + Pj, implies b = 61 or 6 = 62- In this paper we present a novel 
algorithm to compute such atomic fibers. An algorithmic solution to appearing subproblcms, 
computational examples and applications are included as well. 



1. Introduction 

Decomposition of rational polyhcdra is at the heart of several interesting applications. However, 
there are different definitions of decomposability depending on the application. These definitions 
mainly differ in the treatment of the (integer) points of the polyhedron. 

The simplest notion is that of linear decomposition of polyhedra. Two polyhedra P,Q C R™ are 
called homothetic if P = XQ + t for some A > and ( 6 R". Here, a polyhedron P is called 
indecomposable, if any decomposition P = Qi + Q2 implies that both Q\ and Q2 are homothetic to 
P. It can be shown that there are only finitely many indecomposable rational polyhcdra that are 
not homothetic to each other. For further details on this type of decomposition we refer the reader 
for example to Grunbaum (1967), Hcnk et al. (2003), Kannan ct al. (1990), McMullen (1973), 
Meyer (1974), Smilanski (1987). 

Let us now come to a bit more restrictive decomposition. Here we consider only polyhedra of the 
form {1 e R" : Ax < b} for a given matrix A £ Z dxn and varying b G 1 d . To emphasize that we 
only consider integer right-hand sides, we say that a polyhedron P is integrally indecomposable, if 
any decomposition P = Q1+Q2 (into polyhedra with integer right-hand sides) implies that both Q\ 
and Q2 arc homothetic to P. This decomposition is more restrictive than the linear decomposition, 
since only such polyhcdra Q\ and Q2 are allowed that have an integer right-hand side. Hcnk et al. 
(2003) showed finitcness of the system of integrally indecomposable polytopes. This result implies 
important applications: TDI-ness of each member of a family of systems Ax < b, b £ ~L d , can be 
concluded from TDI-ncss of the integrally indecomposable systems. Furthermore the finiteness of 
the system of integrally indecomposable polytopes enables us to compute a finite representation 
of a test set for a mixed-integer linear optimization problem. 

Another important application of integral decomposition of polyhcdra is that of factorizing a 
multivariate polynomial, see for example Abu Salem et al. (2004) and the references therein. Here, 
one considers only polyhcdra of the form {x £ R' 1 : Ax < b} for given matrix A £ Z dxn and 
varying b £ Z d , where each polyhedron is integer, that is, where each polyhedron has only integer 
vertices. Note that the notion of integral decomposability is restricted to integral polyhcdra in this 
application whereas the definition of Hcnk et al. (2003) is valid for arbitrary rational polyhcdra with 
integral right-hand side. The reason for this restriction is the simple observation that the so-called 
Newton polytope Newt(/) := convjai £ supp(/)} associated to a polynomial / = ^2 ieI aix ai 
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with supp(J) = {on : a% 7^ 0} is integer by definition. Moreover, the relation f = gh among three 
polynomials /, g, and h implies Newt(/) = Newt (5) + Newt(h), a theorem due to Ostrowski. 

A direct generalization of the above notion of integral decomposition of integral polyhcdra was 
introduced by Adams et al. (1999). They considered polytopes 

P b := conv{z : Az = b, z £ Z™ } 

called the fibers of b under the linear map /a ■ Z" — > TL d , x 1— > Ax. A fiber Pb is called atomic if 
Pb = + Pb 2 implies b = b\ or b = 62. Note that Pb = P^ + Pb 2 means that every vertex of Pb is 
the sum of a vertex of P bl and a vertex of P b2 (and vice versa). Atomic fibers were used by Adams 
et al. (1999) to construct strong SAGBI bases for subalgebras of polynomial rings. They proved 
that the family of atomic fibers is finite and also gave an algorithm to compute atomic fibers via 
certain standard pairs. Via this algorithm, Adams et al. (1999) computed the atomic fibers of the 
twisted cubic, see Example 2.7. 

In this paper, we consider a slight variation of the notion of atomic fibers that was introduced by 
Maclagan (2001). Instead of considering convex hulls P b of the prcimages 

P b :={z: iz = Me^}, 

of the map /a- Z" — ► Z d , x 1— ► Ax, we consider the prcimages Pb themselves. In Maclagan's more 

general terminology, the sets Pb are called ((0), ^4)-fibcrs; we shall simply call them fibers in the 

remainder of this paper. We call a fiber Pb indecomposable or atomic, if Pb = Pb t + Pb 2 implies 

b = b\ or b = hi- Note that Pb = P bl + Pb 2 means that every lattice point of Pb is the sum of 

a lattice point of Pb t and a lattice point of Pb 2 (and vice versa). This is indeed a very strong 

condition, but again it was shown that there are only finitely many (nonempty) atomic fibers for 

a given matrix A (Maclagan, 2001). Note that atomic fibers are not only minimal (with respect 

to decomposability) within the given family, but also generate every fiber Pb in this family as a 
t 

Minkowski sum Pb = z2 i=1 o^P^, ol% S Z+, where aiP bi stands for iterated Minkowski-addition 
of Pb i with itself. Atomic fibers (of this kind) were used in the computation of minimal vanishing 
sums of roots of unity (Steinbcrger, 2004). 

Recently, the computation of atomic fibers also appeared as a subproblem in the capacitated design 
of telecommunication networks for a given communication demand under survivability conditions 
(Eisenschmidt et al., 2006). In this application, the right-hand side vectors b are taken from 
a sublattice or a submonoid of Z d . Of course, restricting the set of "feasible" right-hand sides 
changes the notion of decomposability, thus the set of atomic fibers is changed. 

Another related notion is that of extended atomic fibers. We call the set 

Q b :={z:Az = b,ze Z"} 

an extended fiber of the linear map of A. We call it atomic, if [Qb fl Oj) = (Qb ± H Oj) + (Qb 2 H Oj) 
holds for all the 2" orthants Oj of R™, then b = bi or b = 62- Here, as well, it can be shown that 
there are only finitely many (nonempty) extended atomic fibers for a given matrix. Also this very 
strong notion of decomposability has an application: the set Woo constructed in Hemmecke and 
Schultz (2003) for use in two-stage stochastic integer programming is in fact the set of extended 
atomic fibers of the family of extended fibers 

{{x, y) : x = b, Tx + Wy = 0, x G T n , y E T 1 } 

where T and W are kept fixed and where b varies. 

Outline. In this paper, we are mainly concerned about designing efficient algorithms for comput- 
ing atomic and extended fibers. The outline of the paper is as follows. In section 2 we first define 
a hierarchy of partially extended fibers that interpolate between fibers and extended fibers. This 
hierarchy not only generalizes the notions of fibers and extended fibers, but also plays a significant 
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role in our algorithms. Motivated by our application in survivable network design, we define de- 
composability with respect to a given finitely generated monoid of feasible right-hand side vectors. 
We prove that, in this more general situation as well, there are only finitely many atomic fibers. 
We also present an algorithmic way to decompose a fiber into a Minkowski sum of indecomposable 
fibers. 

In section 3 we present a first algorithm to compute the atomic extended fibers of a given matrix, 
following the pattern of a completion procedure. We present the algorithm in a simplified setting 
where the right-hand side vectors are restricted to a sublattice (rather than a submonoid) of Z d . By 
restricting the atomic extended fibers to the positive orthants and performing a simple reduction 
step, the atomic fibers (or partially extended fibers) of a matrix can be easily obtained. However, 
this method is not a very efficient one for computing atomic fibers. 

Therefore, we present a more efficient way to compute atomic fibers via a project-and-lift approach 
in section 4 and 5. We present the method in the general setting where a finitely generated monoid 
of right-hand sides is given by its generators. 

Both our algorithms enable us to compute not only the atomic fibers Pb but also the atomic fibers 
Pb according to the definition in Adams et al. (1999). This will be shown at the end of section 3. 

Finally, in section 6, we present first computational results of the project-and-lift algorithm. 

2. (Partially Extended) Atomic Fibers 

Let us now start our treatment with a formal definition of partially extended fibers. 

Definition 2.1. Let AeT dxn be a matrix, b G T d and < k <n. 

(i) The set 

Q { b k) := {z : Az = 6, z G T% x Z"- fc } 

is called an partially extended fiber of order k of the matrix A. The set Qb := Q b * s 
called an extended fiber, and Pb := Qb is called a fiber of the matrix A. 

(ii) Let < I < n. For u,v G R n we say that u C; v i/u w i> w > and \u®\ < \v®\ for all 
components i = 1, . . . , I. We will abbreviate C n by C. For U,V,W C R n we say that 

(i) 

u = v ®w 

and call U the /-restricted Minkowski sum ofV and W, if for all u G U there exist v G V, 
w G W with v,w C ; u and u = v + w. Note that V (B^ W is just the ordinary Minkowski 
sum V + W. We will abbreviate by ®. 

(iii) For < m < n we will denote by n m : R n — > R m with (x\, . . . , x n ) i— » (x\ , . . . , x m ) the 
projection onto the first m components. 

Now we will go on defining atomic partially extended fibers w.r.t. a certain monoid M C ~L d . To 
accompany the hierarchy of partially extended fibers, we define a hierarchy of notions of decom- 
position that interpolates between ordinary Minkowski sums and orthant-wise Minkowski sums. 

Definition 2.2. Let A G Z dxn be a matrix, b G T d and <k,l < n. Additionally, let M C T d be 
a monoid. 

(i) We call Qi atomic w.r.t. and M if there is no decomposition 

with bi,b 2 G M and n{Q$),n{Q$) ^ MQo®)- B V E\ k] (A,M) we denote the set of 
partially extended fibers of order k which are atomic w.r.t. ©W and M. 
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(ii) We denote by E^ (A, M) the set E„ (A, M) and call it the set of partially extended 
atomic fibers w.r.t. the monoid M. We denote by F(A,M) the set E^ n '(A,M) and call it 
the set of atomic fibers w.r.t. M . 

Note that Definition 2.2 also applies to the special case where the monoid M is a lattice. We will 
see later on that it is much easier to compute the atomic (partially extended) fibers of a matrix 
w.r.t. a lattice instead of an arbitrary monoid. 

As our first step, we prove a generalization of the fmitcncss result for the family of atomic fibers. 

(k) 

Lemma 2.3. Let < k < n be fixed. There are only finitely many partially extended fibers Q b 
which are atomic w.r.t. a finitely generated monoid M. 

The proof of this lemma is based on the following nice theorem. 

Theorem 2.4 (Maclagan, 200f). Let k be afield. Let X be an infinite family of monomial ideals 
in a polynomial ring k[x\, . . . ,x n ]. Then there must exist ideals I,JgX with J C J. 

To apply this theorem in our situation of partially extended fibers which are atomic w.r.t. a certain 
finitely generated monoid M , we introduce the following definition. 

Definition 2.5. Let A G Z dxn and M = {mi, . . . ,m t ) C J d a finitely generated monoid. Let 
< k < n be fixed. 

(i) Let a, a G with 

b := ^""^ ^oumi and b := ^ctitrii. 
We say that (a, Qi ) reduces (a, Q b ) and denote 

(a,Q {k) )<l(a,Q ik) ) 

if St C a and Q b = <3g ffi Q b _u- I n particular: b — b G M . 

(ii) We call a pair (a, Q b ) irreducible w.r.t. < if there is no pair (a,Q^) different from 
{a,Q (k) ) and (0,Q^ fe) ) with 

(S,Q {k} ) < (a,Q (k) ). 

Lemma 2.6. Let < k < n be fixed. Let A = {(a 1 , Q^), (a 2 , Qb ), ■ ■ ■} be a set of pairs. 

(i) Let (a\Q b k) ) £ (a^Q (k) ) for all (a*, Q^f), {a j ,Q { b k) ) G A with i < j. Then A is finite. 

(ii) There are only finitely many pairs (a,Q b k ^) which are irreducible w.r.t. <. 

Proof, (i): We associate with a pair (a? , Q b ) the monomial ideal 

n 

I a = (a.(*i.-.*».4+i.*in.-.4»*».^.-.^) : Az = 5^ajmi( = bj), zeZ^x Z n ~ k ) 

i=l 

C Q[ Xl , . . ,,X2n+t], 

where zf = max{0, z{\ and z~ = max{0, — Zi}. Then (a- 7 , Q b ) ^ (a', Q b ) if is not contained 
in I a i . Consider the set X = {I a i, I a 2 , . . .} of ideals associated to the elements in the set A. The 
set X then is an antichain of ideals and is thus finite according to 2.4 (see Maclagan (200f )). The 
finitcness of A follows from the finiteness of X. 

(ii): A pair (a, Q b ) is irreducible w.r.t. < if and only if (a, Q b ) ^ (a, Qr ) for any a ^= a. Let 

A = {(a 1 , Q b k ), (a 2 , Qj* ), ■ ■ ■} be the set of pairs which are irreducible w.r.t. <. Part (i) then 
yields that .A is finite. □ 
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We are now ready to prove Lemma 2.3. 



Proof of Lemma 2.3. It is sufficient to show: for every Q^ "' atomic w.r.t. M there exists a £ J\ 
with 6 = X)i=i a i m i such that (a, Q^) is irreducible w.r.t. <■ Then there is an injective mapping 
from the set of atomic extended fibers into the set of irreducible pairs (a, ) and thus there 
are only finitely many extended atomic fibers w.r.t. M . 

Let b be fixed with an extended atomic fiber w.r.t. M. Let a € Z*|_ with b = Xa=i a i TO i be 
minimal w.r.t. C, i.e., there is no Z*|_ 3a/a with a C a and 6 = &i m i- We claim that the 

pair (a,Q^) is irreducible w.r.t. <. Suppose not. Then there is (a,Q^) < {a,Q^), i.e., a C a 
and = Qg © Q^-l i m Pb/ m g 6 — 6 e M. As is an extended atomic fiber we may w.l.o.g. 
assume that 6 = 6 and 6 — 6 = 0. Therefore 6 = X)'=i Q^ 771 ! an d as a is minimally chosen w.r.t. □ 
we have a = a. This proves our claim. □ 

Example 2.7. In Adams et al. (1999), it was shown how atomic fibers could be used to construct 
strong SAGBI bases for monomial subalgebra over principal ideal domains. As an example, they 
computed the atomic fibers of the matrix -^-=(0123) by hand via an approach different from the 
one we present below. 

In the table below, we list the right-hand sides and all (finitely many) elements in these 18 atomic 
fibers. 
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,(3, 1,1,1), (3, 0,3,0), (2, 2, 2,0), (0,6, 0,0), 



Thus, for example, the fiber given by the right-hand side (8, 7) is not atomic, since it can be 
decomposed into atomic fibers as 

P (?)= P (l) eP (I)- 



This can be quickly verified by looking at the elements in these fibers: 



{(2, 1, 0, 2), (2, 0, 2, 1), (1, 1, 3, 0), (1, 2, 1, 1), (0, 4, 0, 1), (0, 3, 2, 0)} 
= {(0,1,0,1), (0,0,2,0)} ©{(2,0, 0,1), (1,1, 1,0), (0,3, 0,0)}. 
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Indeed, we have 

(2,1,0,2) 
(2,0,2,1) 
(1,1,3,0) 
(1,2,1,1) 
(0,4,0,1) 
(0,3,2,0) 



(0,1,0,1) + 
(0,0,2,0) + 
(0,0,2,0) + 
(0,1,0,1) + 
(0,1,0,1) + 
(0,0,2,0) + 



(2,0,0,1), 
(2,0,0,1), 
(1,1,1,0), 
(1,1,1,0), 
(0,3,0,0), 
(0,3,0,0). 



In Example 2.7 above, it was easy to verify whether a given fiber is a summand in the decomposition 
of another fiber by simply checking the finitely many elements in the fiber for a decomposition. If 
the fibers are not bounded, however, this would not give a finite procedure. The following lemma 
tells us how to solve this problem via the (finitely many!) C-minimal elements in the given fibers. 

Definition 2.8. Let A £ Z dxn and beZ d . Let < k < I < n. 

(k) (k) 

(i) An element v £ Q b is called minimal w.r.t. if there is no w £ Q b with v ^ w and 
w Qi v. 

(ii) We define z,z£ Q b to be equivalent if and only if tti(z) = iri(z). 

For I < n there are infinitely many -minimal elements in general. Therefore we have to restrict 
ourselves to representatives of equivalence classes of Qi-minimal elements. Let R b k ^ denote a 

set of representatives of the equivalence classes of the \—i-minimal elements in Q b . Let these 
representatives be chosen arbitrarily but fixed. 

Remark 2.9. Let A £ Z dx ™ , b £ ~L d and < k,l < n. Then the set of representatives of Q- 

(k) (k) 

minimal elements in Q b , R b , , is finite by the Lemma of Gordan-Dickson ( see for example Cox 
et al. (1992)). 

Lemma 2.10. LetQ<k<l<n and let Q$ ^ 0, Q$ ^ 0. Then Q { b % b2 = © (i) Q$ if and 
only if for every \—i-minimal vector v £ R^ +b2 ; there is a vector w £ Q b ^ with w C ; v. 

Proof. Let v £ Q^ +b . Then there is v £ R b _±_ b ; with v v. Thus, by the assumption in the 
lemma, there is some w £ Q b such that w C ; v C; v. As k < I we have v — w £ ~L k + x Z n_fe and 

/~i(k) - 
thus v — w £ Q bo with v — w C ; v C ; v. 

We now claim that v — (w + v — v) + (v — w) with w + v — v £ Q b , v — w £ Q b , w + v — v Qiv, 
and v — w Qi v, is a desired representation of v. The first two relations are trivial, if we keep in 
mind that Av = Av = b, Aw = b±, b = bi + 62 and k < I. We get the other two relations as follows: 

(a) w + v — v C; v + v — v = v, since by construction tti(w) and tti(v — v) lie in the same 
orthant, and 

(b) v — w v v, since iv C; v. 

Thus, we have constructed for arbitrary v £ Q b ^ +b a valid representation of v as a sum of two 

elements from Qj® and Q b whose projection onto the first I components lie in the same orthant 
as the projection of v onto its first / components. This concludes the proof. □ 



Using this lemma repeatedly, we are now able to find, for a given right-hand side b £ M, a 
decomposition Q b = ©j =1 &iQ b k \ cxi £ Z+, that is, we can find a decomposition of a partially 
extended fiber into a sum of partially extended fibers which are atomic w.r.t. the monoid M . 
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Algorithm 2.1 Algorithm to decompose extended fibers into sums of extended atomic fibers 

Input: A, right-hand sides {by, . . . , b s } of the set of extended atomic fibers E^ k \A, M) 

Output: ai, . . . , a a such that Q b = (J) aiQi 

1: oi\ := . . . := a s := 

2: for i = 1 to s do 

3: while Q ( b k) = © Ql%. and b - b t G M do 

4: b := b-bi 

5: a 4 := a, + 1 

6: end while 

7: end for 

8: return: a-y,..., a s . 



It remains to state an algorithm that computes the finitely many C-minimal elements in Q b for 
fixed k. We will do this in the following paragraphs. 

We have to find for some I G {1, ...,n} and some k G all C; -minimal elements in 

(projections of) fibers of the form 

MQ { b k) ) = {(x,v) eZ k + x T {l - k) : 3z G T {n - l) with A(x,y,z) = b}. 

If b = 0, then is the only C ; -minimal element. If not, we reduce this problem to the problem 
of finding a Hilbert basis of a cone. It is not hard to show that all Cj-minimal elements (x,y,z) 
correspond to the elements (x, y + , y~, z, 1) in a Hilbert basis of the cone 

{(0, y+,y-, z, u) G Z"+ ( '- fc ) +1 : A(x, y+ - y~,z) - bu = 0, x, y+, y~,u> 0}. 

In general, this is not a pointed rational polyhedral cone (and thus need not have a unique inclusion- 
minimal Hilbert basis), since there can be linear relations among the (free) variables z. However, 
projected onto the space of the variables x, y + , y~ , u, the nonnegativity constraints lead to a pointed 
rational polyhedral cone that possesses a unique inclusion-minimal Hilbert basis. Such a minimal 
Hilbert basis can be computed for example with 4ti2 (see 4ti2 team). 

Note that the splitting of y into y + and y~ is only used for exposition here. In practice, one can 
directly use y when computing the IZ;-minimal elements, see Hemmecke (2006, Section 2.6) for 
more details. 

3. COMPUTATION OF (EXTENDED) ATOMIC FIBERS 

In the following we show how to compute the finitely many (extended) atomic fibers of a matrix 
A G Z mx ™ w.r.t. a lattice A. In this section we will present a simple algorithm; we will give a 
more complex and much more efficient algorithm in the following sections. Both algorithms use 
the algorithmic pattern of a completion procedure. 

We will denote the columns of matrix A by Ax,...,A n G Z m . Note that the function 
normal form(s,G) in Algorithm 3.1 stems from Algorithm 3.2. 

Lemma 3.1. Algorithm 3.1 terminates and computes a set G such that {Q^ & ■ b G G} contains 
all atomic fibers of A. 

Proof. Associate with b G A the monomial ideal I am '■= {x^ z+ ' z ** : Az = b 7 z G Z n ) C 
Q[x\, . . . , X2n], where (z + ) J = max(0, z°) and {z~)° = max(0, — z°) for all components j = 1, . . . ,n. 
Algorithm 3.1 generates a sequence {/i, / 2 , . . .} in G\F such that Qf j ^ Q/ ; © Qfr-fi whenever 
i < j. Thus, the corresponding sequence {IaJi , Iaj 2 , • • •} of monomial ideals satisfies Iaj, ^ IaJi 
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Algorithm 3.1 Algorithm to compute extended atomic fibers 
Input: F := {±61, . . . , ±b s } with {b x , . . . , b s ) = A n AT 1 

Output: a set G, such that {Qb '■ b £ G} contains all extended fibers of A which are atomic 
w.r.t. A 
1: G := F 

2: C := |J {/ + g} (forming S-vectors) 

3: while C ^ do 

4: s := an clement in C 

5: C:=C\{s} 

6: / := normal form (s, G) 

7: if / ^ then 

8: G:=GU{/} 

9: G := G U Ugeci/ + 3} (adding S-vectors) 

10: end if 

11: end while 

12: G:=GU{0} 

13: return: G. 



Algorithm 3.2 Normal form algorithm 
Input: s, G 

Output: a normal form of s with respect to G 
l: while there is some g € G such that Q s = Q g © Q s _ g do 
2: s := s — g 
3: end while 
4: return: s 



whenever i < j. We conclude, by Theorem 2.4 given by Maclagan (2001), that this sequence of 
monomial ideals must be finite and thus, Algorithm 3.1 must terminate. 

It remains to prove correctness. For this, let G denote the set that is returned by Algorithm 3.1. 
Moreover, let Qi be an extended atomic fiber of A with 6^0. We will show that b £ G. 

Since F \ {0} C G \ {0}, we know that Qi = J^Qbj f° r finitely many (not necessarily distinct) 
bj e G \ {0}. This implies in particular, that every z £ Qi can be written as a sum z = with 
Vj £ Qbj ■ We will show that we can find vectors bj £ G such that every z £ Qi can be written as a 
sum z = v j with Uj £ Q^. and «j C z. This implies Qi — Q Qbj ■ Since Qi is atomic, and thus 
indecomposable, this representation must be trivial, that is, it has to be = Q^, and therefore 
we conclude b £ G. 

With Lemma 2.10 it is sufficient to consider the C-minimal elements in Qi, = {zi, . . . , z^}, to 
decide if it decomposes w.r.t. ©. From all representations Qi = ^2j £ j Qb., with bj £ G\{0} choose 
a representation and elements Vij £ Qbj with Zi = Ylj^j v i,j * = 1, . . • , fe, such that the sum 

»=i jeJ 

is minimal. By the triangle inequality we have that 

EEM^EiMii- (2) 

i=i j'e j i=i 
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Herein, equality holds if and only if all Vij have the same sign pattern as Zi, i = 1, . . . , fc, that 
is, if and only if we have Vij C z^ for all i and all j. Thus, if we have equality in (2) for such a 
minimal representation Q b = Yljcj Qbj, then Vij G Qbj and Vi_j C Zi for all occurring Vij, and we 
are done. 

(It should be noted that we have required bj G G\ {0} for all appearing bj, that is in particular, 
bj 0. Those bj will be sufficient to generate all C-minimal elements in the extended fiber Qb- 
We get the remaining elements in Qi by adding elements from Q$.) 

Therefore, let us assume that 

EE Mi >£ INI*- ( 3 ) 

i=i jej i=i 

In the following we construct a new representation Qi = Ylj'ej'Qb'. an( i elements v\j whose 
corresponding sum (1) is smaller than the minimally chosen sum. This contradiction proves that 
we have indeed equality in (2) and our claim is proved. 

From (3) we conclude that there are indices io, ji,]2 an d a component m G {1, . . . ,n} such that 



v 



(m) (to) 



l 0,Jl 10,32 



v) { < 0. As bj, , bj, G G, the sum 6,, + 6,-, has been built and the extended fiber Q bi +&, 



has either been reduced to Qo by sets Qb , ( , j" £ J", during the Algorithm 3.2 or bj 1 + bj 2 has been 
added to G. In the latter case we set J" := {j"} with by := + bj 2 . This gives representations 



+ u «,j2 = "^i.j" with ,j" G Qbj,, and i3 -» C v. 

for i = 1, . . . , fc. As all Wi.j» lie in the same orthant of R n as Vi_j 1 + Wi j2 , we get 



E 

j"e,7" 

with strict inequality for i = ig. 



= IKji + ^,j 2 Hi < IKjilli + lk J2 ||i, 

l 



Thus, replacing in Qb = J2jejQbj the term Q bj± + Q bj2 by J2j»t=j"Qbj„, w e arrive at a new 
representation Qi = Ylj'eJ' Qbj, whose corresponding sum (1) is at most 

k k 

EE Ki'iii< EE Mi> 

i=l j'GJ' i=l jeJ 

contradicting the minimality of the representation Qi = Ylj^ j Qbj- This concludes the proof. □ 

Remark 3.2. One may of course use Algorithm 3.1 also for the problem of finding indecomposable 
extended fibers among the elements in the family of extended fibers Qb = {z : Az = b, z G Z™} 
where A is kept fixed and where b is allowed to vary on the lattice which is spanned by the columns 
of matrix A. The input set then becomes F = {±^4i, . . . , ±A n }. 



Having an algorithm available that computes all extended atomic fibers w.r.t. a given lattice A, we 
can use it to compute partially extended atomic fibers w.r.t. © and A: If is atomic then so 
is Qb, as any decomposition of Qb, restricted to x Z n ~ fc , would give a decomposition of Q[ . 
This way of computing partially extended atomic fibers of a given matrix A G Z mxra is illustrated 
in Figure 1 and formalized in Algorithm 3.3. 

The solid arrow from the bottom up in Figure 1 stands for the completion procedure which is 
given by Algorithm 3.1. The dashed arrow from the top to the bottom illustrates the procedure 
of intersecting the extended atomic fibers with Z5. x J_ n ~ k and dropping the reducible (or empty) 
fibers afterwards. 
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Algorithm 3.3 Computing partially extended atomic fibers 
Input: F := {±61, . . . , ±b s } with (b x , . . . , b s ) = A n AT 1 , keT+ 

Output: A set G* such that {Q^ : b S G*} contains all partially extended fibers of order k 
which are atomic w.r.t. © and A 
l: Apply Algorithm 3.1 to the set F. Let G denote the output. 
2: G* := 0. 

3: for b e G with Qf ] ^ do 

4: if Q ( b k) ^ Q { g k) © Qf} g for all g ± b e G then 

5: G* := G* U {b} 

6: end if 

7: end for 

8: return G* 





{Q b } 


b£G 


completion / 




restriction, 

dropping 

reducible 


{QbjbeF 


1 

{Qt^beG* 


+ 


© 



Figure 1 . Computing (partially extended) atomic fibers via extended atomic fibers 

Being given the atomic fibers Pb of a matrix A it is easy to compute the atomic fibers Pb which 
have been defined in Adams et al. (1999). Recall that Pb := conv{z : Az = b z G Z™} and that 
Pb is said to be atomic if each decomposition P = Pb x + Pb 2 implies either b = b\ or b = 62. 

Lemma 3.3. If Pb is an atomic fiber of the matrix A then Pb is atomic, too. 

Proof. Suppose Pb = P^ +Pb 2 ( an d &1, ^2 7^ 0). Then we have: Pb = conv(Pj,) = conv(P(, 1 +Pb 2 ) = 
conv(Pf >1 ) + conv(P;, 2 ) = Pf >1 + Pb 2 which is a contradiction. □ 

Lemma 3.3 enables us to compute the atomic fibers Pb via Algorithm 3.4. 

4. Preliminaries of the project-and-lift algorithm 

The way of computing atomic fibers presented in section 3, however, is pretty slow, since there 
are far more extended atomic fibers than atomic fibers. A similar behavior can be observed when 
one extracts the Hilbert basis of the cone {x : Ax = Q,x £ R™} from the Graver basis of A, as 
the Graver basis is usually much bigger than the Hilbert basis one is interested in. Hcmmcckc 
(2002) showed that one can reduce this difference in sizes by a project-and-lift algorithm. With 
this algorithm, bigger Hilbert bases, even with more than 500,000 elements, can be computed 
nowadays. 
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Algorithm 3.4 Computing the atomic fibers Pb 

Input: F := {61, . . . , b s } with Pb i is an atomic fiber 

Output: A set G = {bi, . . . , bt] such that P\. is an atomic fiber. 

l: Set G := 0. 

2: for all b £ F do 

3: if P b ^P g + P b -g for all b ^ g £ F then 
4: G:=GU{b} 
5: end if 
6: end for 
7: return G 



In this section and in the following one, we will present a similar algorithm to compute the atomic 
fibers of a given matrix A £ Z dxn which is significantly faster than Algorithm 3.3. This algorithm 
puts us in the position to compute not only the atomic fibers of a matrix but the atomic fibers 
w.r.t. an arbitrary (finitely generated monoid) M, i.e., the right-hand side b is only allowed to 
vary in this monoid. During the algorithm we consider partially extended fibers ^ = {z £ 
Z'jT 1 x z™- fc +i . _ £,j w ith varying b £ M w.r.t. /c-restricted Minkowski-sums. 

Let M = (mi, . . . , m t ) C Z d be a finitely generated monoid and let A £ ~Z dxn be a matrix. We 
want to compute the atomic fibers of matrix A w.r.t. the monoid M. The algorithm proceeds in 
n individual steps. The fc-th step is illustrated in Figure 2. 



k-1 


k 




\Q b HeGk-i 
1 1 

1 {Qf }b£F k 


T / 




(k-1) 


(fc) 
© 



Figure 2. The fc-th step of the project-and-lift algorithm 

The fc-th lifting step follows the arrows in the figure. It starts by performing a "preprocessing step" 
in which the input set is prepared for the main part of this lifting step. This process is illustrated 
by the dotted arrow and will be explained in more detail in section 5.3. 

The fc-th lifting step continues as follows: it performs a completion step similar to the one we 
presented in Algorithm 3.1, which is illustrated by the solid arrow going from the bottom up. This 
step will be explained in more detail in section 5.1. 
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The dashed arrow, finally, stands for a step where we drop all elements of the fibers having a 
negative k-th component. It might happen that an atomic partially extended fiber becomes empty 
or reducible when processing this last step. Therefore we have to perform another reducibility test. 
The details of this subroutine will be given in section 5.2. 

Having performed the fc-th lifting step we continue performing the (k+ l)-st lifting step. The whole 
project-and-lift algorithm is illustrated in Figure 3. After having performed n of these lifting steps 
we arrive at the finitely many fibers of the matrix A which are atomic w.r.t. M. 



\k 





1 


2 


3 




n 






{Qb° } }beG 




{QT}bEG 2 
1 1 
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/; 

{Qb htt. 


/] 
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/ 1 
/ 1 

/ ! 

/ {QbheG* 
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Wfc >bGF 


T / 

{Qb {bePt 


T / 

W b Sb£F 2 
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W b Sb£F 3 








(0) 


(1) 
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(2) 
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(3) 
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(n) 
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Figure 3. The scheme of the project-and-lift algorithm 



Dealing with infinitely many atomic fibers. Let A G J dxn be a matrix and M C Z d a monoid 
which is finitely generated. The project-and-lift algorithm will deal with partially extended fibers 
w.r.t. ©W and the monoid M where I < n. Recall from Definition 2.2 that Qff is atomic w.r.t. ©^ 
and M if there is no decomposition 

Qi k) = Off § Off 

with fei , b% G M and 717 (Qff ) , 717 (Qff ) 7^ 717 (Qq ) . Note that for I < n there are usually some 
b G M with 7r;(Qff ) = TTiiQ^). Therefore, if Q b k] is atomic w.r.t. ©W and M then so is Q^, 
^1+26' ' ' ' ' means that for I < n we usually have infinitely many partially extended fibers 

which are atomic w.r.t. ©*■''. It is clear that no terminating algorithm may compute the whole 
set of atomic partially extended fibers w.r.t. ©^ and M. Therefore we introduce a preorder <i 
(i.e., a reflexive and transitive binary relation) on the set of right-hand side vectors b G M with 
non-empty partially extended fiber and perform the l-th step of the project-and-lift algorithm 
w.r.t. the preorder 

Definition 4.1. Let C M be the submonoid of M with Q ( b k) ^ for be M (fc) . Let A G Z dxn , 
< k < I < n and let b,b G M^. We say that b±bifb-b€ where = {\ l+1 A l+1 + 
... + X n A n : A, 6 Z}fl M . 

b G M ^ is called minimal w.r.t. if there is no b^be AfW with b ^1 b. 



COMPUTATION OF ATOMIC FIBERS OF Z-LINEAR MAPS 



13 



Note that b^ib implies m(Qf ') = ni(Q b ) and Q b > = Qf> ©W Q b \- 

The relation -<i defines a preorder on the set of right-hand sides b £ M with non-empty partially 
extended fibers of order fc. Additionally we have the following relation between the sets §"': 

{0} = S* (n) C S (n - 1] C...C§° = M {0) . (4) 

Lemma 4.2. Let M ~ {mi, . . . , m t ) be a monoid which is finitely generated. 

(i) Let < k < I < n and let M D F = {b\, bi, . . •} be a set of vectors with bi bj for all 
i < j . Then F is finite. 

(ii) Let < fc < I < n and let F — {b\, bi, ■ ■ ■} be a set of right hand sides satisfying is 
atomic w.r.t. and M and bi ^1 bj for all bi ^ bj. Then F is finite. 

Proof, (i): Let bi,bj G F with i < j and let 0^,0? € with bi = X)l-=i a k m k and bj = 
Sfe=i Q; i TO 'i:- Then (a^Q^) ^ (a- 7 , Q^). Suppose not. Then we have a 1 C a J and Qj? fc ' = 
Qb\ ®Qb^-bi wn i cn implies that Q b ' = Qb\ ®^ Qb -b,- But this last relation contradicts the 
fact that bi fa bj. Therefore {a\Q b ^ ) ^ (a j ,Q { b f) for all b u bj € F for i < j. Finiteness of F 
follows with Lemma 2.6 (i). 

(ii): This is a direct consequence of (i). □ 

Our algorithm will work with sets of vectors F which have the property claimed in Lemma 4.2. 
Additionally they will admit the following property: if b £ M is the right-hand side of a partially 
extended fiber which is atomic w.r.t. ©^ and M then there is b G F with b ^1 b. This means 
in particular: lib E M is minimal w.r.t. and is atomic w.r.t. ©^ and M then 6 e F. Note, 
however, that the converse is not true in general: It is not guaranteed that for every b 6 M there 
is a b ^1 b that is minimal w.r.t. 

5. The k-th step of the project-and-lift algorithm 

In the following subsections we will explain the individual steps the project-and-lift algorithm 
performs during one lifting step. 

5.1. The completion procedure. In this subsection we will explain the so-called "completion 
procedure" in the fc-th step of the algorithm. This part is illustrated in Figure 4. 

Let M be a monoid which is finitely gener ated. We denote by = {b e M : Q { b k) ^ 0} the 

submonoid of all right-hand sides b S M having non-empty partially extended fibers of order fc. 
We have: 

Af D Af (0) DlW D ... D M {n \ 
As in the previous section, , < I < n, will denote the set {Aj+i^4;+i + . . . + X n A n : \i £ Z}flM. 

Definition 5.1. We introduce a weight function uj m for partially extended fibers Q b (m < fc < n) 
by 

Urn{Q ( b k) ) =min{||7r m (u)||i: v G Q^}. 

Remark 5.2. Actually, it suffices to determine ||7T m (u)||i for \— m -minimal elements v in Q b to 
determine the value of UJ m (Q b ). To see this, suppose there is w S non-minimal w.r.t. C m . 

Then there is v G Q b with v \— m w and thus < v 3 < w 3 for j = 1, . . . ,fc. Therefore |-7r m (u)||i < 
I km (w) ||i. 

Lemma 5.3. Algorithm 5.1 with input set F^-i := {61, . . . ,b s } and monoid terminates 
and computes a set Gk-i = G" =0 U G"- 1 U {0} C M^ k ~^ with properties (i) and (ii). 
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Figure 4. The completion procedure of the fc-th lifting step 
For the proof of Lemma 5.3, we have to introduce some more notation. 

Notation 5.4. During the proof of Algorithm 5.1 we examine elements of partially extended fibers. 
These elements will be denoted as follows: 

Q^~ l) 3 z = (21, Z2, *a) G Z^ 1 x Z x Z"~ fc , 

i.e., z\ £ 1 denotes the first k — 1 components, Z2 £ Z i/ie /c-i/i component and z% £ lf n ~ k ^ 
denotes the last n — k components. 

We will use the following lemma in the proof of Lemma 5.3. 

Lemma 5.5. Let F fc -i C M^" 1 ) foe a sef admitting the following property: for every right-hand 
side b £ of a partially extended fiber Q^ k 1 ' which is atomic w.r.t. ©( fc ~ 1 -' and M there 

exists b £ Fk-i with b <k b. Let [3 £ M^ -1 ' be the right-hand side of an arbitrary partially 
extended fiber of order (k — 1). Then we find bi £ F^—i such that for M^^ 1 ^ 3 f3 := ^2 bi we have 
dik P and 

«r"-er«T"- <*> 

Proof. Let [3 £ AI^^ 1 ^ 1 \ {0} be the right-hand side of a partially extended fiber of order (k — 
1). Consider a decomposition of into a sum of fibers of order (k — 1) which are atomic 

w.r.t. ©( fe_1 ) and M: 

(fe-l) 

Q { r ] - e Qt i] 

As the partially extended fibers ^ are atomic w.r.t. ffi( fc_1 ) and M there are 6, £ Ffc_i with 
6, ^ k b, for all i £ I. Consider M f*" 1 ) 3 /3 := £. t£/ 6 4 . We have /3 ^ fc /3 because fo* - £ £ 
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Algorithm 5.1 The completion procedure to compute atomic partially extended fibers w.r.t. a 
monoid 

Input: A set Ffc_i C A/^ -1 ) with the following properties: 

(i) For every right-hand side b G M^"' -1 ) of a partially extended fiber Q^ -1 ^ which is 
atomic w.r.t. Q^ -1 ^ and M there exists b G Fk-i with b rk- b. 

(ii) bi bj for 6*, &j G F fc _i with bj ^ b, y 
Output: A set Gk-i Q M^ k ~ x ' with the properties: 

(i) For every right-hand side b G M^" -1 ) of a partially extended fiber which is 
atomic w.r.t. ®^ and M there exists b G Gk-i with & b. 

(ii) &j ^ fc 6j for bi ^ ^ G G k -i. 

l: G-=° := {/ G F fc _! : ^(Q^^) = 0} 

2: C"=° :=U/ >9e G"=o{/ + ff} 

3: while G w =° ^ do 

4: s := an element in C* w=0 

5: G w=0 := G w=0 \ {s} 

6: / := monoid-normal-form (s, G w=0 , 0, M( fc_1) ) 

7: if / i then 

8; G w=0 := G w=0 U {/} 

9: G"=°:=G"=°U|J ge( >=o {/ + <?} 
10: end if 

11: end while 

12: G w=0 := 

13: for all b G G u=0 do 

14: if Q { b k - 1] ^ Qf~ l) 8 (fe) Q^ x) for all b ^ g G G w=0 with b — g £ M then 
15: G w=0 := G w=0 U {&} 
16: end if 

17: end for 

18: G^ 1 := {/ G F fc _x : w fc _i(Q^ _1) ) > 0}, G^ 1 = 

19: for all g G G^ 1 do 

20: G"^ 1 = G^U monoid-normal-form (g, G"=°, 0, M^ -1 )) 

21: end for 

22: G"- 1 :=U/, aeG ">i{/ + fl} 

23: while G"^ 1 ^ do 

24: s := an element in G"- 1 with smallest weight u>k-i(Qs k ) 

25: G"^ 1 := G"^ 1 \ {.s} 

26: / := monoid-normal-form (s, G w=0 , G<^\ M^" 1 )) 

27: if / ^ then 

28: G*^ 1 := G*^ 1 U {/} 

29: G"- 1 := G"- 1 U U fl6G — ug<^ {/ + .9} 

30: end if 

31: end while 

32: G k -i ■= G w=0 U G^ 1 U {0} 

33: return Gk-i 



implies J^iei h-bi = (3 - f3 e S^. Additionally 
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Algorithm 5.2 The monoid-normal-form algorithm 

Input: s, G u=0 , G"- 1 , membership oracle for M (fe_1) 
Output: a normal form of s w.r.t. G w=0 U G"- 1 and M^ -1 ) 

l: if 3g G G"^ 1 with Qi fc_1) = Q [ g~ 1] @ {k) and s-ge M^" 1 ) then 

2: return 

3: else 

4: while 3g G G w=0 with Qi fe_1) = Q ( g k ~ 1] ©(*) Q^~g ] and s - g e M^- 1 ) do 
5: s := s — g 
6: end while 
7: return s 
8: end if 



which together with /3 = J2iei ^ implies that Q ~ fe ^ = Q~ fc and our claim is proved. □ 
We are now in the position to proof Lemma 5.3. 

Proof of Lemma 5. 3. As the following proof will be slightly complex, consider the following outline 
of the proof first. 

(1) We show that G fc _i C M^" 1 ). 

(2) We show that the set G" =0 is finite and that for all h,bj 6 G u=0 with 6, ^ 6 we have 

(3) We show that Gk-i is finite. This implies that Algorithm 5.1 terminates. At the same 
time we show that the output set admits property (ii), i.e., bi bj for bi bj G Gk-i- 

(4) We show that if is an atomic partially extended fiber w.r.t. and M then there 
is b <k b with b G Gk-x- This is property (i) of the output set. 

Step 1. 

It is clear that Algorithm 5.1 returns a set Gk-i C M k ~ l . This is guaranteed by the monoid- 
normal- form algorithm, where we ensure that the elements added lie in M^ k ~ 1 ' . 

Step 2. 

We will now prove that the set G" =0 is finite. To this aim we show finiteness of G w=0 first. 
Consider the sequence G w=0 \ {/ G Ffc-i : Uk-iiQ^ 1 ^) = 0} = {fx, f%, . . .} produced in lines 
1-11 of the algorithm. Clearly fi G M^ k ~^ for all i. Additionally fi -£k fj for a U i < j- Suppose 
not and let f t < k fj. Then Q {k ~ 1] = © (fc) Qf.Z^- As fj has been added to G w=0 , the 

second criterion of the monoid-normal-form algorithm is not satisfied, i.e., fj — fi £ M*-' 0-1 -'. But 
fj — fi S S^ implies in particular that fj — fi G M and as Qj.^f. ^ we have fj — fi G M( fe_1 ) 
which is a contradiction. Therefore fi ^ fj for all fi, fj G {/i,/2, ■ ■ ■} with i < j. Finiteness of 
G w=0 follows with Lemma 4.2. 

As G"=° C G" =0 it is clear now that G w=0 is finite. Additionally lines 12-17 of Algorithm 5.1 
guarantee that bi ^fe bj for all 6j, 6j G G w=0 with 6j 7^ bj. 

Step 3. 

Let G W=Q := {b G G fc _i : 0J k -i{Q { ^ 1] ) = a} for a G Z+. Furthermore let G^ Q := {6 G G fc _i : 
^k-iiQ^ ) < a} for a G Z+. We will show via induction that that bi ^-k bj for G G"- a 
with 6j 7^ bj. Lemma 4.2 then yields that G"- a is finite. Clearly we have Gk-i = U Q ez + G- a . 
Let G Gfc_i with bi =/= bj. Then there is a G Z+ with bi,bj G G"^ Q yielding b t ^ k bj. The 
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set Gk-i admits property (ii) of the output set thus which together with Lemma 4.2 yields that 
Gfc_i is finite. 

We will show via induction that G"- a is finite. With step 2 of the proof we know that our claim 
is proved for a = 0. Suppose that our assertions are true for all integers lower or equal than a. 
We will prove our claim for a + 1. Let bi,bj £ Q^<a+i an( j SU pp 0SC \ >i -< fc There are several 
cases possible: 

(i) bi,bj £ F k _x 

This contradicts input property (ii) of the input set F^-i- 

(ii) h £ F fc _i, bj £ Fa.-i 

This contradicts the if-clause of Algorithm 5.2 because bi then is an appropriate reducer 
of bj._ 

(hi) h <£F k -i, bj £ F fe _i 

As Wk-i(Qi k Zl] ) = and as G" =0 is completed before G"- 1 we know that there is b £ G w=0 
with b -<k bj — bi <k bj. But this is a contradiction to lines 18-21 of Algorithm 5.1 because 
in this case bj would not have been added to G"- 1 then. 

(iv) bi,bj £ Ffc-i 

Depending on whether either bi has been added to G W = Q + 1 befre bj was added or not we 
cither have a contradiction to the if-clause of Algorithm 5.2 or to the else-clause of this 
algorithm. 

We know via induction that G"- Q admits property (ii) of the output set. We will now show that 
this is also true for G u '~ a+1 . Let bi, bj £ G"- a+1 and suppose bi <k bj. By induction, the previous 
discussion and monotonicity of the weight-function Wfc(-): 6j £ G w - Q and bj £ G W=Q+1 . But this 
contradicts the if-clause of Algorithm 5.2. Therefore Q UJ < a + 1 admits property (ii) of the output 
set which had to be proved. 

Step 4- 

Let b £ M( fc_1 ) such that Q^* is atomic with respect to (B^ and M. With Lemma 5.5 we know 
that there is b £ M, b <k b, admitting a representation (5): 

where bi £ F^-i- We will show that b £ G^-i- The above representation implies in particular that 
every z = (zi,z 2 ,z 3 ) £ 1) can be written as (zi,z 2 ,z 3 ) = Yliei i z h z \i z l) wit; h ( z h z h z l) e 
Q^ k . In particular: 

{zi,4,4)^k-l(zuZ2,Z 3 ) (7) 

for all i. If Q^ k 3 (z\,z\,zX) C fe (zi,Z2,z 3 ) was valid this then would imply: ^ = 

m (fc) n (k-i) 

Let i?j 1 ' = {(z\, z\, 4)i • ■ ■ j z 2i ^3)} ^ e the set of representatives of the □^-minimal elements 
in Q~ fe ^ according to Definition 2.8. With Lemma 2.10 we know that it is sufficient to analyze 
the IZ fe -minimal elements in a partially extended fiber to decide decomposability w.r.t. @^- k \ 

From all representations Q^~ 1 ^ = ©ygj Qj fc— ^ with bj £ Gk-i and ^k(Qf ) k ~ 1 " > ) =/= 
n k (Q^ ^) and where the C fe _ 1 -minimal elements in R^ k ^ are represented as (z\,z l 2 ,z^) = 
Sje/ ( z i J : Z 2 J ) Z 3 J ) w ith ( z i J i Z 2 J ) z 3 ^ ) S G;^ 1 for i = !,...,£, choose a representation and 
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elements (z 1 , z 2 , z 3 ) such that the sum 



EEii(^^2 J )iii (8) 

1=1 jeJ 



is minimal. 



By the triangle inequality we have 

EEii(^^2 J ')ii^Eii(^'4)iii ( 9 ) 

i=i je.J i=i 

Herein equality holds if and only if all {z{ 3 , z 2 ' 3 ) have the same sign pattern as (z\, z 2 ), i = 1, . . . , t, 
that is if and only if we have (z^ 3 , z 2 3 , z 3 3 ) (z{,z 2 ,z 3 ) for all j G J and all i = 1, 
Thus if we have equality in (9) for such a minimal representation Q~ fc l ' = ©jgj 1 ^ ^ then 
by Lemma 2.10 Qi fe_1) = ©^Q^ 1 ' and as 7r fc (Q^ _X) ) ^ 7r fc (Q£ fc_1) ) and as Q[ fe_1) and thus 



j(fe-i) 
'S 

Therefore let us assume, that 



Q~ is atomic w.r.t. ©( fc ) and M this representation must be trivial and we are done. 



/ 



EEiiK J ' z 2 J )iii>Eii( z "i^2)iii (io) 

i=l jeJ i=l 

In the following, we construct a new representation Q~ fe ^ = j? Qj? and elements 

(,S 1 '' 7 , ^2 J > ^ ) whose corresponding sum (8) is smaller than the minimally chosen sum. This con- 
tradiction proves that we have indeed equality in (9) and our claim is proved. 

From (10) and from (7), i.e., z\' 3 C z\ for all i and j, we conclude that there are indices 
such that z^ ' 31 ■ z l 2 '' n < 0. As bj ± , bj 2 G G u=0 the sum bj 1 +bj 2 was built during the algorithm. We 
have u>k-i(Qfy~+l ) = ^fe-iCQ^ ) + w fc-l(Q&^ ^) — an d thus there is no partially extended 
fiber with weight Wk-i greater than which reduces Q^ k _^ . Consequently the partially extended 
fiber was either reduced to Q fc by sets , j" G J", during the else-clause of the 

monoid-normal- form algorithm or the vector bj 1 + bj 2 has been added to the set G w=0 . Then either 
b n +b j2 G G w=0 or we find sets Qf~ X) ,j" G J", with - ©J-k/" with V e G< " =0 - 

In the former case, set J" := {j"} with 6j« := bj 1 + bj 2 . 

This gives representations 

/ 2,71 i.jl 2,71 \ , / 2.72 2,72 2,72\ \ ^ /~2,7 

(z x J1 , z 2 ' , z 3 ' ) + ( Zl J - , z 2 ' , z 3 ' ) = 2^ ( z i 

j"eJ" 

(^l'"' , S '•' , Sg'"' ) G Q5 .„ ^ , (^l'"' 1 2 2 '^ , Za' 3 ) Efc (^1 J , Zn/ 3 , Z3'"' ) + {Z{ 3 , Zj'"' 7 2 '3'"' ) 

for i = 1, . . . ,t. As all {z\ 3 , Zj' 7 ) he in the same orthant as {z\' 31 , z 2 31 ) + (z\' 32 , z 2 32 ) we get: 



,z 2 ,z 3 



E ) 

i"eJ" 



1 ^2 



with strict inequality for i = £q. Thus, by replacing in 1 ' = ©jgj 1 '' Q^" 1 the term 
Q (k-i) ffiCfc _ 1} g(fc-i) by 0(^ (( g(fc-i) wc arrivc at a new represcntation q^- 1 ) = 0<*-$ Q^ 1} 

whose corresponding sum (8) is at most 

EE ii(^'^F)iii<EEii(^'^)iii 

■i=i j'e.T 2=1 jeJ 
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contradicting the minimality of the representation Q~ fe 1 ' = ©jgj Q^. ■ Therefore we have 
equality in (9) and thus b G G k -i concluding our proof. □ 

5.2. Intersecting with the appropriate orthant and testing reducibility. In this subsec- 
tion we want to illustrate the step of the project-and-lift algorithm which follows the completion 
procedure in each lifting step. This "intersection and reducibility test" is illustrated by the dashed 
arrow in Figure 5. 



k 




(*■) 



Figure 5. Intersecting with the appropriate orthant and dropping reducible par- 
tially extended fibers 



Algorithm 5.3 Intersecting and testing reducibility 
Input: A set Gk-i C M*- fc_1 ) with the properties: 

(i) For every right-hand side b G M^ k ~^ of a partially extended fiber Q b k which is 
atomic w.r.t. ©( fc ) and M there exists b G G k -i with b < k b. 

(ii) bi bj for bi,bj G G k -i with b { ^ bj 

Output: A set F k C of right-hand sides with: 

(i) For every right-hand side b G of a partially extended fiber Q b which is atomic 
w.r.t. and M there exists b € F k with b ^ k b. 

(ii) bi bj for bi, bj G F k with 6, ^ bj. 
1: F k := 

2: for all b G Gfc_i with Q { b k) ^ do 

3: if Q ( b k) ^ Q { g k) © (fc > Q { b } g for all 6 ^ .g G G fc _i with b-g&M then 
4: F fc := ^ fe U {6} 
5: end if 
6: end for 
7: return F k 



Lemma 5.6. Algorithm 5.3 with input set G k -\ terminates and computes a set F k C with 
the properties (i) and (ii). 
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Proof. Termination of Algorithm 5.3 is clear. We have to show correctness of the algorithm. 
But this is easy as well: if 6 G M^ 5 " 1 ' and Q (k) ^ then b G AfW. Therefore F k C Af«. 
If Qtf* is atomic w.r.t. ®^ and M, then Q^ k ^ is atomic w.r.t. (B^ and M as well, because 

Qb°^ Q Qb" an d every decomposition of Q^ k ^ would give a decomposition of . This 
characteristic immediately implies property (i) of the output set because we have property (i) of 
the input set. 

To see property (ii) of the output set, suppose that there are 61,62 G F^ with 62 61. Then, 
= ®( fe ) Q[^L &2 and 61 - 6 2 G §W . In particular, 61 - 6 2 £ A/ which is a contradiction 
as 61 would not have been added to Fk in this case. This yields that bi ^ k bj for all 6,,6j G -Ffe. 
Therefore Algorithm 5.3 is correct and terminates. □ 

5.3. Refining the preorder. There is one more step to explain in the fc-th lifting step of the 
project-and-lift algorithm. This step is illustrated by the dotted arrow in Figure 6; it is implemented 
in Algorithm 5.4. 



k 


{Qfh 




\ 


r 


(fe) 

© 



Figure 6. Refining the preorder to prepare the k + 1-st lifting step 
Lemma 5.7. Algorithm 5.4 terminates and is correct. 

Proof. Termination of the above algorithm is clear once we have shown that we can construct a 
finite set C with the property that 

Vs g S (k) 3s {j) g C with s w r<fc+i s. 

For this, let us first construct generators for the monoid 5'' fc+1 - ) . These can be found by considering 
the homogeneous system of linear equations 

n t 
j=k+2 r=l 

in the variables s G Z d , a G Z+, and in A G Z" _fc_1 . If we extract the values of s for all finitely 
many minimal homogeneous solutions of this linear system, we obtain a generating set {si, . . . , s p } 
for the monoid §( k+1 > . 
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Algorithm 5.4 Refining the preorder 

Input: A set Fk C Afw of right-hand sides with: 

(i) For every right-hand side b G of a partially extended fiber Qj® which is atomic 
w.r.t. and M there exists b G F k with b b. 

(ii) b t ^ fc bj for ^,6, G F k 

Output: A set Fk C Af( fc ) with the following properties: 

(i) For every right-hand side b G Af' fc ) of a partially extended fiber which is atomic 
w.r.t. © (fc ) and M there exists 6 G F k with 6 ^ fc+ i b. 

(ii) 6j ^fe+i 6j for all bi, bj G Fk with 6, ^ 6j. 
l: Compute a set £ = {s«, . . . ,s^} C S« with: 

Vs G S w 3s® G £ with ^fc+i s 

2: SctF fe := U U {& + *}• 
3: for all & G F k do 

4: if 36 G Fk with 6^6 and b d?k+i b then 
5: F k := Fk \ {b} 
6: end if 
7: end for 
8: return Fk 



Now let us consider the finite set 

F = XiS,: < Xi < l,i = n {Afe+iTlfc+i + . . . + X n A n : A l G Z}. 

For each / G F we now consider the set (/ + 5( fc+1 )) CiS^ and construct a finite set £/ of vectors 
in (/ + 5 (fe+1) ) n such that 

Vs G (f + 5 (fc+1) ) n 5 (fc) 3s (j) G C f with s w ^fe+i s. 

Then £ = U/gf — ^ s finite and has the desired property. 

In order to construct £/, let us consider the inhomogeneous system of linear equations and in- 
equalities 

n t n t 

f + XjAj = f + a r m r = s = A; = (3 r m r 

j=k+2 r=l i=fc+l r=l 

in the variables s G Z rf , a, (3 G Z+, and in A G Z n ~ fc ~ 1 ,/i G Z"~ fe . The left-hand part states 
s G / + 5( fc+1 ) and the right-hand part encodes s G S^ k \ Then a suitable set Cf can be found 
by computing the finitely many minimal inhomogeneous solutions to this linear system and by 
collecting the corresponding values of s. We have thus proved that we may construct a finite set 
£ admitting the claimed property It remains to prove that the set Fk constructed from £ admits 
the properties claimed in Algorithm 5.4. 

It is clear that F k C M W because F k C AfW and £ C 5^ C M. Now let 6 G A/ (fc) be the 
right-hand side of an atomic partially extended fiber w.r.t. ®' fc ) and M. Because of property (i) 
of the input-set, we find b G Fk with b < k b, i.e., there is s G S^ k ' with b = b + s. As s G there 
is s G £ with s r^fc+i s implying that there is s G S^ k+1 ^ with s = s + s. Let b := b + s. Then 
clearly 6 ^fc+i 6. Thus either b G or there is b' G Fk with &' ^^+1 6 dfc+i b. This concludes the 
proof. □ 



COMPUTATION OF ATOMIC FIBERS OF Z-LINEAR MAPS 



22 



5.4. Initial input set and final output set. The previous subsections have shown how one 
step of the project-and-lift algorithm works. We have to perform n of these steps to obtain a set 
G* which contains the right-hand sides of all atomic fibers w.r.t. the monoid M = (mi, . . . ,m t ). 
We start with the initial input set Fq = {0}. This is a valid input set because every non-empty 
extended fiber with right-hand side b G M, i.e., b G M^°\ also lies in S^. 

5.5. Simplifications for the lattice case. As already mentioned at the beginning of section 4, 
the project-and-lift algorithm to compute the atomic fibers of a matrix with right-hand side b 
varying on a lattice A is much easier to treat than the case of general monoids. 

The simplifications of the project-and-lift algorithm are based on the fact that the difference of two 
arbitrary lattice vectors &i, hi G A, b\ — b 2 , is again a lattice vector. This fact has implications for 
the preorder -<i on the right-hand side vectors b G A. Let A G Z dx ", A C ~i_ d a lattice and < I < n. 
Consider the preorder introduced in section 4: bx b 2 if b 2 — bx £ §"> = {Xi + \Ai + i + . . . + X n A n : 
Xi G Z}nA. As A is a lattice we obtain that 6i —62 € S"> as well. This means: b\ -<i b 2 b 2 d?i ki- 
ln other words, the preorder is in fact an equivalence relation. 

Our aim in this subsection is to simplify the refining step in our project-and-lift algorithm. Recall 
from section 5.1 that the input set F^-i satisfies the following two properties which ensure finiteness 
and correctness of the algorithm: 

(i) For every right-hand side b G M^ -1 ' of a partially extended fiber Q^ k which is atomic 
w.r.t. ©( fc_1 ) and M there exists b G F k -i with b < k b. 

(ii) bi bj for b { ^ bj G G k -x- 

In this subsection, we will define a new input set F^-x of Algorithm 5.1 which may be computed 
much easier than the set F k _\. Having defined this new input set we will expose some properties 
of it. Finally we will show that the new input set Fj^—x is sufficient to guarantee finiteness and 
correctness of the completion procedure, i.e., of Algorithm 5.1. 

Let Ff~-x be the output set of Algorithm 5.3 and consider the following integer program: 

min A 

k 

S.t. X k A k + = E Pjh 

i>k+X j=X 

Xk-X < (11) 
-A fc -A < 
A > 1 

A, Aj, yUj G Z 

There are two possible cases: either the integer program (11) is infeasible or it admits an optimal 
solution X*,\*,n1j. Consider the former case first and let bx, 62 & A with bx dk-x b 2 . Then b\ < k b 2 . 
This is the case because b 2 — bx G A and b 2 — 61 = J2i>k ^« 6 ^. If the absolute value of 
Afc was greater or equal than 1 the difference b 2 — b\ would imply a feasible solution of the integer 
program (11). Thus X k = and therefore we have b\ < k b 2 . In this case we set F k _\ := F^-x- 

Now consider the latter case. Let s := Y^i= k K^i- Again there are two possible cases: either 
s d?k —s or s ^ k — s. In the former case we set Fk-x := F k -i U {s}, in the latter case we set 
F fc _x := F fe _i U {±s}. 

Lemma 5.8. Let (11) admit an optimal solution A*,A*,/i*. We assume w.l.o.g. that A£ G Z+. 
Let be atomic w.r.t. ©( fc_1 ) and A. Then there is b G -F&— 1 and At, G Z with b + Xb-s d k b. 

Proof. Let fee A with Q^ 1 ' atomic w.r.t. and A. Then there is b G F fc _i with b ^ fc _i 6, 

a consequence of Fk-x being the output set of Algorithm 5.3. Let s := b — b = J2i>k ^i^i- Then 
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there is n G Z and < f < A£ with A& = n - \* k +f. We will show that r = 0. To this aim consider 
s := s — n ■ s = fAk + J2i>k+i (^» — n ' K)Ai- As s, s G A we have so e A. If f ^ the feasible 
solution of (11) implied by sq admits an objective value lower than A*, because f < A£. This 
contradicts the optimality of the solution A*, A*, fi*. Therefore f = and thus b + n ■ s b. □ 

Lemma 5.9. For bi, bj G -Ffc-i with bi ^ bj we have bi -£k bj- 
Proof. There are a few different cases to consider: 

Case 1: bi,bj G Fk-i- Then bi -Ak bj as F^-i is the output set of Algorithm 5.3. 

Case 2: bi G F/-— i and bj = s. We have to show that s bj. Suppose not and consider 
bj — s = E t >Hi ^^'»- But tncn bj = Ei>it+i + V! , • ' ; an d thus ^fc_i which 
contradicts bj being an element of Fk—i. 

Case 3: bi = s and bj = — s. Here we have s y^k ~ s by our assumptions. □ 

Lemma 5.9 implies that F^—i defined as above satisfies property (ii) of the input set of Algorithm 
5.1. We continue giving another property of the set Fk-i- 

Lemma 5.10. Let b G A k ~ 1 . Then there is b -<k b with 

(fc-i) 

Qf ~ l) = © Qf" 15 where b t G F k -l- (12) 

i 

Proof. As b G A (fc_1) there is b d?k-i b with 

(fc-i) 

Ql~ 1) =®QZ~ 1) where bj G Fk—i- 

j 

This representation is a consequence of Ff-—i being the output set of Algorithm 5.3. If the integer 
program (11) is infcasiblc then b b and our claim is proved. Therefore let (11) admit an optimal 
solution. W.l.o.g. we assume that F^-i = Fk-\ U {±s}. Consider a decomposition of Q^ 1 ^ 
into a restricted Minkowski sum of partially extended fibers which are atomic w.r.t. this restricted 
Minkowski sum and the lattice A: 

Q { b ~ 1] = © QbT^- (13) 

i 

With Lemma 5.8 we know that for each bi there is 6, G -Ffe-i and Ab i G Z such that bi + X^ -s ^ bi. 
We set b := J2i bi + • s. Then we have b <kb and 

This proves our claim. □ 

Now we want to show that the input set Fk-i is sufficient to guarantee finiteness and correctness 
of Algorithm 5.1. An input set admitting properties (i) and (ii) is sufficient to do so. We have seen 
in Lemma 5.9 that our set Fk-i admits property (ii). It admits property (i) as well if the integer 
program (11) is infcasiblc. But it does not admit this property in general if the integer program 
(11) is feasible. Note that in the proof of Algorithm 5.1 property (i) is only used to guarantee 
a representation (5) according to Lemma 5.5 with the projection of the summands satisfying 
Kk(Qb ) Tk(Qo )• But with Lemma 5.10 this representation may be guaranteed as well. 
Furthermore i^kiQ^ ) ^ i"k(Qo )■ This is clear for b G Fk-i because ^fc-i b. As ^ s we 
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Algorithm 5.5 Refining the preorder (equivalence relation) in the lattice case 
Input: A lattice A = (li,.-.,lt) and a set Fk Q A^ of right-hand sides with: 

(i) For every right-hand side b £ A^ of a partially extended fiber which is atomic 
w.r.t. 0( fe ) and A there exists b £ F k with b ^ k b. 

(ii) bi bj for b it bj £ F k with b t ^ bj 

Output: The set Fk Q A^ defines as above. 
1: Solve the following integer program: 

min A 

k 

S.t. \ k+ iA k+ i + = E 

i>k+2 j=l 

A fe+ i-A < (14) 
-A fe+ i-A < 
A > 1 



A, Aj,yUj £ Z 

if (11) is feasible then 

Let A*, A*,/i* be an optimal solution of (11). 
Set s := E,> fc+ i KU 
if s <k+i —s then 

return F k := F k U {s} 
else 

return F k := F k U {±s} 
end if 
else 

return F k := F k 
end if 



have an analogue result for n k (Q s ). This finally implies that Algorithm 5.1 terminates and is 
correct when given input set F k -\. 

Besides the modification of the input set, Algorithm 5.1 stays the same. Of course we may drop 
all tests if b — g £ A during the normal-form algorithm because for b, g £ A it is clear that the 
difference b — g £ A. The same is valid for Algorithm 5.3. It stays the same except for the dropping 
of tests whether b — g £ A. 

Algorithm 6 is substituted by the above Algorithm 5.5. It does not compute the set F k but the 
set F k . 

Lemma 5.11. Algorithm 5.5 terminates and is correct. 



Proof. This is a direct consequence of the discussion in this subsection. 



□ 
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6. First computational results 

We have created an implementation of the "project-and-lift" algorithm for the lattice case (sec- 
tion 5.5). The implementation is written in Allegro Common Lisp 8.0 and C. For the computation 
of the minimal elements of partially extended fibers, we use the library libzsolve, which is a part 
of 4ti2 (4ti2 team), version 1.3.1. In this section, we report on the computational experience with 
this code on several test problems. All computation times are given in CPU seconds on a Sun Fire 
V440 with UltraSPARC-IIIi processors running at 1.6 GHz. 

6.1. Results for number-partitioning problems. We first consider the problem of partitioning 
a natural number n into given parts (natural numbers) a\, . . . , a k (with possible multiplicity). To 
this end, consider the set 

k 

P n = { (xi, . . . ,x k ) e : n = ^Xi • Oj j. (15) 

We are interested in a minimal set {n-i, . . . , n q } of natural numbers such that the set P ra of partitions 
of every number n is the Minkowski sum of some of the sets P nj . Thus we are interested in the 
atomic fibers corresponding to the matrix 

( ai a 2 a 3 • • • a k ) . (16) 

We consider this problem for various sets of numbers a±,. . . ,a k - The results are shown in Table 1. 

Table 1. Results for number-partitioning problems 



Parts 


Atomic fibers 


Time (s) 


1 


1 


1 


1 2 


2 


1 


1 2 3 


4 


1 


12 3 4 


9 


1 


1 2 3 4 5 


32 


875 


12 3 4 5 6 


41 


>1000 


2 3 


3 


1 


2 3 5 


14 


1 


2 3 5 7 


72 


149661 


3 5 


1 


1 


3 5 7 


30 


1 



6.2. Results for homogeneous number-partitioning problems. Next we consider the prob- 
lem of partitioning a natural number n into given natural numbers a\ , . . . , a k (with possible 
multiplicity), where we prescribe the number of summands. To this end, we consider the set 

k k 

p n = { (xi,...,x k ) el+in^^srfl,, to = J2 X * }' 

i=l i=l 

We are interested in a minimal set {(rai,ni), . . . , (m q ,n q )} of pairs (m,n) such that the set P™ 
of partitions of every number n into to summands is the Minkowski sum of some of the sets P™/ ■ 
Thus we are interested in the atomic fibers corresponding to the matrix 

/i i - i 

\ ai a 2 ■ ■ ■ a k 

Again we consider the problem for various sets of numbers oi,. . . The results are shown in 
Table 2. We remark that the problem data (1,2,3,4) correspond to a problem equivalent to the 
one from Example 2.7. 
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Table 2. Results for homogeneous number-partitioning problems 



Parts 


Atomic fibers 


Time (s) 


i 

X 


l 

X 


1 

X 


1 9 


9 

Z 


1 

X 


1 9 % 

X Z O 


\ 

-i 


1 

X 


1 9 1 A 

X Z O ^ 


1 X 
1 O 


1 

X 


X Z O T: (J 


7Q 


1 Q 

X J 


1 1 Q H 

x z o o 


1 9 


1 

X 


1 9 9 fi 
X A o u 


o o 


9 

Z 


1917 

X Z O 1 


1 Q 


1 
X 


X A o o 


■JO 


•-to 


1 9 1 Q 

1 Z O J 




o 
z 


X Z O XU 


o / 


zuu 


X Z O X X 


t.) ,j 


u 


x z o x z 


^ 99 


xuzu 


1 2 3 13 


52 


21 


1 2 3 14 


163 


5136 


1 2 3 15 


67 


72 


1 2 3 17 


79 


216 


2 3 


2 


1 


2 3 5 


4 


1 


2 3 5 7 


26 


1 


2 3 5 7 11 


262 


152792 



6.3. Results for Steinberger's sums of roots of unity. One example that appears and was 
solved in Stcinbcrgcr (2004) is the computation of the atomic fibers of the matrix 



/ 1 


-i 





-1 


1 











°\ 





i 


-1 





-1 


1 




















1 


-1 





-1 


1 





V o 











1 


-1 





-1 


1 / 



This matrix corresponds to a certain problem on 3 x 3 tables and has in fact 31 atomic fibers and 
79 extended atomic fibers. The atomic fibers can be computed with our implementation in less 
than one CPU second. 

The next higher problem on 4 x 4 tables leads to the matrix 



/ 1 


-1 








-1 


1 





























M 





1 


-1 








-1 


1 



































1 


-1 








-1 


1 






































1 


-1 








-1 


1 



































1 


-1 








-1 


1 



































1 


-1 








-1 


1 






































1 


-1 








-1 


1 



































1 


-1 








-1 


1 





1 





























1 


-1 








-1 


1 J 



Our implementation was able to compute the 12, 675 atomic fibers for this matrix within 6.5 CPU 
days on a Sun Fire V440 with UltraSPARC-IIIi processors running at 1.6 GHz. 
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